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ABSTRACT 



We present a statistical model which shows the influence of turbulence on a thermonuclear 
flame propagating in C+0 white dwarf matter. Based on a Monte Carlo description of 
turbulence, it provides a method for investigating the physics in the so-called distributed burning 
regime. Using this method we perform numerical simulations of turbulent flames and show that 
in this particular regime the flamclct model for the turbulent flame velocity loses its validity. In 
fact, at high turbulent intensities burning in the distributed regime can lead to a deceleration 
of the turbulent flame and thus induces a competing process to turbulent effects that cause 
a higher flame speed. It is also shown that in dense C+0 matter turbulent heat transport is 
described adequately by the Peclet number, Pe, rather than by the Reynolds number, which 
means that flame propagation is decoupled from small-scale turbulence. Finally, at the onset of 
our results we argue that the available turbulent energy in an exploding C+0 white dwarf is 
probably too low in order to make a deflagration to detonation transition possible. 



Subject headings: methods: statistical - nuclear reactions - stars: supernovae: general - 
turbulence 
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Introduction 



The thermonuclear explosion of a C+0 Chandrasekhar-mass white dwarf, which is believed to be the 
underlying process of a type la supernova (SN la), has been subject to numerous investigations. However, 



despite the fact that there are different plausible models explaining the history of the explosion (Nomoto 



et al. 1984; Woosley fc Weaver 1986; Miiller fc Arnett 1986; Livne 1993; Khokhlov 1991a; Arnett and 



Livnc 1994; Niemeyer fc Hillebrandt 1995; Hoflich, Khokhlov fc Wheeler 1995; Hoflich 1995; Wheeler et 



al. 1995 ; Hoflich fc Khokhlov 1996 ), many important details still remain unclear. Thermonuclear reactions 
provide the source of energy which possibly unbinds the white dwarf. Thus, a SN la is characterized by the 
physics of thermonuclear flames that propagate through the star. The physical conditions of this flame vary 
drastically during the different temporal and spatial stages of this process. A better understanding of these 
conditions, their interaction with the flame and finally their consequences regarding the explosion itself are 
still important issues from the theoretical point of view. 

Here, we focus on the interaction between thermonuclear burning and turbulence taking place during 
the burning process. Turbulence is caused by different kinds of instabilities (like shear instabilities or the 



Rayleigh- Taylor instability) that occur on certain length and time scales (for an overview, see Niemeyer 
fc Woosley 1997). Rough estimates give a turbulent Reynolds number Re « 10^^ at an integral scale of 



L w 10^ cm ( Hillebrandt fc Niemeyer 1997 ). Consequently, the Kolmogorov scale Zk, i-e- the scale where 
microscopic dissipation becomes important, is about 10^* cm. This enormous dynamical range makes a 
direct representation - at least by means of numerical methods - practically impossible. On the other hand, 
turbulence is a characteristic feature during the explosion process and it must be considered in any realistic 
model of a SN la. The simultaneous coupling to energy generation due to nuclear reactions leads directly 
to the physics of turbulent combustion, where a few crucial problems, even for terrestrial conditions, still 
remain unsolved. These problems include the prescription of the effective turbulent ffame speed or the 
existence of different modes in turbulent combustion along with their physical properties. Global features 
of turbulent combustion can be systematically classified by a small set of dimensionless parameters. This 
classification, which already found a wide utilization among the chemical combustion community, can be 
used in the field of SNe la in order to emphasize universality in the physics of turbulent combustion. 

In this work we study the properties of turbulent flames in a certain state, the so-called distributed 
flame regime. The physical characterization of the latter is given by the situation where turbulent motions 
are fast enough to disturb the flame on microscopic scales. Therefore a distributed flame does locally not 
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look like a laminar flame anymore which is a basic difference to turbulent flames in the flamelet regime. 
Turbulent burning fronts in an exploding white dwarf are flamelets at high and intermediate densities, but 



below ^ 5 X lO^g cm one cannot expect that the flamelet picture is still valid ( Nicmeyer fc Woosley 1997 ). 

Until now, in the context of SNe la, burning in the distributed regime has only been considered 
qualitatively. Therefore we wish to present first quantitative results regarding this issue. Furthermore, 
there are good reasons to believe that this combustion mode plays an important role in the deflagration to 
detonation transition (DDT), which in turn is a promising model, for empirical and theoretical reasons, for 
the explosive stage of a SN la (e.g. Khokhlov 1991h| ; Nicmeyer fc Woosley 1997 ; Nicmeyer fc Kerstein 1997| ; 



Khokhlov et al. 19971) . 



In order to attempt a representation turbulent dynamics in the distributed flame regime, we use a 
new model, formulated in one spatial dimension, which nevertheless provides essential features of three 
dimensional homogeneous turbulence. It consists of a statistical description of turbulent mixing and a 
deterministic evolution of the underlying microphysics. This method allows a systematical investigation 
of turbulence phenomena. When it is coupled to a nuclear reaction network, it gives first insights of how 
the flame structure is affected by turbulence on scales, which have not been resolved in direct numerical 
simulations. In particular, we investigate the flame properties in cases where the Gibson scale is comparable 
to the thickness of an undisturbed conductive flame. The Gibson scale is defined as the length scale 
on which the turbulent velocity fluctuations equal the laminar flame velocity. In the case of flames in 
degenerate white dwarf matter Iq becomes comparable to the thickness of the flame only for densities 



around 2 x 10^ g cm ^ and below ( Khokhlov et al. 1997 ; Nicmeyer fc Woosley 1997 ). 

The outline of the paper is as follows: first we introduce and describe the statistical method that we 
use to model turbulence. Then we couple the latter to physically relevant microscopical diffusion processes, 
such as temperature and viscous diffusion, and to external energy sources coming from nuclear reactions. 
The resulting method is eventually used to investigate some properties of turbulent flames in the distributed 
regime, like their effective flame speed. Finally we discuss the relevance of our results to the DDT problem. 



2. One-Dimensional Turbulence 



Since fundamental aspects of turbulence can be recovered from the knowledge of the statistical 
moments and correlations of the velocity flow, the statistical approach to turbulence is particulary 
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appealing. Therefore, the history of statistical methods in turbulence theory is rather long. In this context 



we present a novel model of turbulence (Kerstcin 199£). It is a stochastic method, realized as a Monte 
Carlo simulation, which allows to compute statistical properties of the flow velocity and of passive scalars 
in stationary and decaying homogeneous turbulence. One-dimcnsional turbulence (ODT) represents many 
aspects of three-dimensional turbulence, but it is formulated in only one spatial dimension. This model 
provides the temporal evolution of a characteristic transverse velocity profile u{y,t) of the turbulent 
medium, where y is the spatial location on a finite domain [0, Y] and t is the elapsed time. This is done 
in a two-fold way: u{y,t) is subject to a molecular diffusion process and to a random sequence of profile 
rearrangements representing turbulent eddies. Reflecting the typical behavior of turbulence kinematics, the 
event rate of these profile rearrangements (so-called eddy mappings) is proportional to a locally averaged 
shear of the velocity profile u. 

Given the transverse velocity profile at a certain time, u{y,t), the mapping which models the action of an 
turbulent eddy of size I and at the position yg reads 

"(3y - 2yo, t) yo<y <yo + fil 

uiy,t) = ■{ui-3y + 4yo + 2l,t) yo + fil < y < yo + f2l (1) 
u{3y~2yo~2l,t) yo + f2l < y < yo + I 

where it is /i = 1 — /2 = Vs and u{y,t) = u{y,t) for all y ^ [yo,z/o + This three-valued map represents 
the typical features of a turbulent vortex, namely rotation and compression. Each random map defines an 
eddy time scale, T{yo,l,t), via 

2 A\ui{yo,t) - ui{yo + l/2,t)\ 

Here, ui is the boxcar averaged profile of u over a length scale I. A is the only model independent and 
dimensionless parameter which has to be fixed empirically. Vortical kinetic energy is fed by the kinetic 
energy of the local shear. Thus equation can be interpreted as an energy balance. It then is 

i Po {I T-^f = i Po (2 A \ui{y^, t) - uiiyo + 1/2, t)\)\ (3) 

where po is the fluid density. Equation (|^) resp. (|^) is used to introduce the statistical hypothesis of ODT. 
It assumes that the occurrence of eddies with size I and location yo (with respective tolerances dl, dy) is 
governed by a Poisson random process with mean event rate 

—-^-—dydl^:X{yo,l,t)dydl, (4) 
l'^T[yo,l,t) 
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where the processes for different values of yo and I are statistically independent. In this context A can be 
viewed as a factor which scales the event rate A. Motivated by the results in Kerstein (1999) we choose 
A = 0.23. In addition, the microscopic evolution is modelled by a diffusion process of the kind ut = i^Uyy, 
with the kinematic viscosity v. The temporal evolution of turbulence, stationary or decaying, depends on 
the choice of boundary conditions for u. A decaying turbulent intensity is obtained by periodic boundary 
conditions, whereas the stationary case is given by the choice of jump periodic boundary conditions: 
u{y + Y,t) — u{y) + uq, uq ■— u{Y, t ^ 0) and u{y, t = 0) being strictly monotonic on [0, Y]. 

The numerical implementation of ODT reproduces many typical features of three-dimensional 



homogeneous turbulence ( Kerstein 1999 ). For instance, power spectra of the kinetic energy show the self 
similar k~^/^ power law within the inertial range down to the scale where the transition to dissipation takes 
place, c.f. Fig. 1. An obvious advantage of this ansatz is the high spatial resolution of turbulence compared 
with multidimensional numerical models. In combination with the relatively moderate computational effort, 
ODT appears as a useful tool for performing parameter studies in turbulence theory. On the other hand, 
ODT does not consider any pressure fluctuations (dynamical or external) in the temporal evolution of the 
velocity profile u. The reason for this artefact is the inherent conservation of kinetic energy in only one 
spatial velocity component due to equation (^. Note that the pressure gradient term in the Navier-Stokes 
equations redistributes energy among the different spatial components. No such redistribution is considered 
here. However, for isobaric flows ODT appears to be an appropriate model of turbulence. 



3. Turbulent Flames in Dense C+O Matter 

Although the aim of this work is to investigate some properties of turbulent burning in dense C+O 
matter, we begin this section by recalling some properties of laminar, conductive flames. Since the physics 
of undisturbed flames inside a white dwarf is relatively well understood, the first reason for it is to verify 
established results, such as those given in the work of Timmes & Woosley (1992). The other reason is that 
thermal conduction and viscous diffusion will be employed as the underlying microscopical processes for 
ODT in order to model turbulent burning fronts in white dwarfs. 

The state of unburned matter, i.e. density and nuclear composition, uniquely defines the propagation 
velocity s\ of the conductive flame. Timmes & Woosley (1992) calculated the flame velocities as well as their 
thickness for various fuel compositions and densities. At lower densities, for lO^gcm^'^ ^ P lO^gcm"^, 
the speed of the laminar flame decreases rapidly. This behavior is accompanied by a strong increase of 
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the flame thickness S\, which is essentiafly the size of the nuclear reactive zone (Timmes & Woosley 1992 



Khokhlov ct al. 1997| ). As aheady mentioned in the introduction, we are especiaUy interested in densities 



of the order ~ (1 — 3) x lO^gcm^^. Thus, we set up a conductive flame propagating into unburned matter 
consisting of half ^^C and half ^^O at densities of p = 1.3 x 10''gcm^'^ and of p = 2.3 x 10''gcm^^. This 
is done by solving the equations for the conservation of mole fractions and enthalpy in planar geometry, viz. 

dY 

-^^Y. -^'^kXjk (i) + YiYkXkj (i) , (5) 

j-k 

— - —— (a—^ - P—- + —S (6) 

dt pcp dx \ dx J dt p Cp ' 

i 

Herein S denotes the local specific energy generation rate, a the thermal conductivity and Na is Avogadro 
number and Bi is the nuclear binding energy of the nucleus considered. The nuclear reaction network 
consists of seven species, viz. '^He, ^^C, ^^O, ^"Ne, ^^Mg, ^'^Si, and ^^Ni. Equation (||) already makes use of 
the fact that the Lewis number, that is the ratio of heat diffusion to mass diffusion, in white dwarf matter is 
around 10^ and consequently microscopic transport of the element species can be neglected. Furthermore, 
in white dwarf matter the laminar ffame speed is much smaller than the speed of sound, si < O.OOlag. 
Therefore pressure does not change significantly across the ffame front. 



Finally, the typical time scale for free collapse or free expansion of the star can be estimated by (Fowler 



fc Hoyle 1964) 



446 

Th-^s. (8) 



P 

For a white dwarf of radius i?wd = 10^ cm and of mass A/^d = 1.4 Mq, equation gives « 0.02 s. 
Timmes & Woosley (1992) stated that gravitational infiuence can be safely dropped since the diffusion 
timescale of the ffame, « 5\ / D, where 5\ is the ffame thickness and D is the temperature diffusion 
coefficient, is several orders of magnitude smaller than Ta_. However, for very low densities (p <^ lO^gcm""^) 
these two timescales can become comparable and a possible expansion of the star could directly affect the 
laminar ffame. To avoid additional complexity we assume in this work that the hydrodynamical timescales 
is much longer than the relevant timescales of the burning front. Thus we consider the small-scale burning 
front to be not affected by stellar expansion. 

The numerical solution of equations (^) and (|^) is given explicitly in time combined with an appropriate 
equation of state representing ions, black body radiation, electrons and positrons. For p = 2.3 x lO^gcm^'^, 
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Figure 2 shows the propagation of the laminar flame which is indicated by the moving jump in the 
temperature profile. In the steady state evolution, the flame moves into the unburned matter at a velocity 
s\ ~ 2.1 X fO^cm s^^. The flame thickness Si is characterized by the width of the temperature jump 
caused by carbon burning. It is 5i = 1.1 cm. After ^^C destruction the burned material has a temperature 
3.2 X 10^ K and a density of 1.45 x 10^ g cm"'^. For these conditions the subsequent oxygen burning has a 
destruction timescale of Tdcs,o ~ 0-1 s. That means that as the carbon flame propagates with a steady-state 
velocity sj, oxygen burning ignites at a distance of 2.1 x 10"^ cm behind the flame front. This distance is 
about ten times larger than our whole computational domain. Consequently, we are not able to track the 
whole spatial region where nuclear reactions take place that are triggered by carbon burning. The extension 
of the reactive zone is even larger if one takes burning of heavier elements, such as silicon, into account. 
Thus, the resolution of the flame with all its reactions up to nuclear statistical equilibrium would require 
an immense spatial resolution. 

However, in our case the fuel consists half of carbon and the destruction of the latter contributes nearly 
all of the total energy release by the laminar flame. In addition, the speed of the laminar flame is governed 
by carbon destruction. Therefore it is fair to concentrate only on the latter. We also solved the above 
equations for a density of 1.3 x 10^ g cm^'^ with the same nuclear composition as before. The resulting 
flame velocity for this case is s\ = 0.85 x lO^cm s~^, while the flame's thickness is 3.8 cm. Our results agree 
well with the ones given in the work of Timmes & Woosley (1992). 

As already mentioned in the introduction, the Gibson scale Iq can be used to measure the influence 
of turbulence on a flame. Since it is directly related to the turbulent fluctuation velocity u' , the Gibson 
scale depends strongly on the kind of turbulence that is considered. Essential features of turbulence, like 
structure functions, the temporal evolution of the turbulent intensity and the energy cascade, depend 
on the absence or presence of certain geometrical (homogeneity, isotropy) and physical (external forces, 
energy sources, state of the turbulent matter) conditions. On large scales in the interior of a white dwarf, 
turbulence is mainly caused by the Rayleigh- Taylor instability of the flame which takes place on a length 
scale of /rt ~ 10® cm. This instability produces a turbulent energy cascade down to dissipative length scales 



Ik « (10""* - 10^3) cm dKhokhlov 1995| ; [Niemeyer fc Woosley 19971 ). On aU scales Ik < I <^ Iu.t turbulence 
is believed to be decoupled from gravitational effects and thus can be described by the Kolmogorov theory. 
The latter sets the scaling law for u' to be u'{l) ^ Inbetween these scale bounds nuclear burning is 



affected by turbulent motion and different kinds of burning regimes reveal (For an overview, see Niemeyer 



fe Woosley 1997). Here, we put our emphasis on the special situation where homogeneous, isotropic 
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turbulence interacts with the microscopic structure of the laminar flame, i.e. where Iq becomes comparable 
to ^1. The relation Iq « S\ marks the transition into the distributed burning regime, where the smallest 
turbulent vortices can enter the interior of the laminar flame and can carry away reactive material before 



it is completely burned ( [Peters 1986| ; [Niemeyer fc Woosley 1997] ). Taking into account the small Prandtl 



numbers in the dense matter of a white dwarf, we can describe this situation by a constraint on the 
turbulent Karlovitz number. This number is the ratio of the diffusion timescale and the eddy-turnover-time 
at the smallest length scale of turbulent heat transport. The latter turns out to be the Kolmogorov scale, 
Zk, only in case of Pr « 1. For Prandtl numbers much smaller than one turbulent heat transport is subject 
to significant diffusion effects already at a scale of ZkPr Thus we have 

Ka.^PrV^^l. (9) 
Using Kolmogorov theory an equivalent representation of the latter condition reads 



u'{L) Sl 

ri 

Re-i/2pr-i/2>i. (10) 



Ka = PrV^Re^/^ 
Sl L 



There is a simple relation between the Gibson length and the Karlovitz number as follows 

Ic^SiKa-^. (11) 

Therefore, Iq < 5\ is equivalent to Ka > 1. The actual value Ka for a given turbulent flame plays an 
outstanding role in turbulent combustion physics. 

We now combine ODT with the structure of a laminar flame to show how certain properties of a 
nuclear flame in dense degenerate matter of white dwarf depend on the Karlovitz number. The obvious 
step is to incorporate equation (^ together with a viscous transport evolution for the transversal velocity 
profile u{x, t). Then turbulent advection is modelled by the random eddy mappings generated by ODT, c.f. 
equation (0). This system reads as 

Ou_P. f^dp_\ ^^^^ 



dt pcp dx \ dx 



— = 1 d (dT\ ^dl/p ^ 1 ^ 
dt pc-p dx \ dx j dt Cp ' 



where m, p, T are the rearranged profiles of the velocity, the density and the temperature according to the 
eddy- mapping, c.f. equation (^. 
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We choose the amphtude of velocity fluctuations of u to be ~ 10^ cm s^^ at ~ 10^ cm, and to 
obey Kolmogorov scahng. These values are motivated through the expected speed of buoyant unstable 
hot bubbles of size L 10^ cm, which become Ray leigh- Taylor unstable and eventually give the main 
contribution to turbulent energy on large scales. Having fixed this velocity at a certain lengthscale, we 
can use Kolmogorov scaling to estimate u'il) on smaller scales I. Thus the amplitude of turbulent velocity 
fluctuations in our model is comparable to the expected small-scale velocity fluctuations within a SN la. 

However, due to computational limitations we cannot consider u to be the actual fluid velocity of the 
white dwarf matter. Since its Prandtl number is around Pr J> 10^^, the numerical resolution of both, the 
viscous and the heat conductive diffusion scale, would cost too much numerical effort even in one spatial 
dimension. To get out of this dilemma, we use an 'artiflcial' Prandtl of 2 x 10^^ and of 10^^. This choice 
underestimates the inertial range of turbulent motions within an exploding Chandrasekhar-mass white 
dwarf. But we believe that the physics of distributed burning in this kind of matter is decoupled from 
scales much smaller than ZkPr~^^^, because turbulent eddies of sizes well below this limit will be smeared 
instantly out by temperature diffusion. Therefore, we do not essentially change the flame properties in the 
distributed burning regime by increasing Pr. Of course, this is only valid for sub-unity Prandtl numbers, 
otherwise one possibly leaves the distributed regime. 

The random velocity fleld u is generated by stationary ODT, thus u' becomes constant in time after 
an initial transient phase. For our studies we perform simulations where u' (at a length scale / — 1 cm) is 
equal or larger than a given value s\. Then we estimate some flame properties such as the effective turbulent 
flame speed, st- So far, analytical results do not predict the hmction ST(u';Ka) correctly. For instance, 
experimentally observed bending effects or even complete flame quenching are still mostly unexplained 
problems (Ronney 1995). 

An overview of the various numerical models is shown in Table |^. Each realization is uniquely 
parametrized by its Reynolds, Karlovitz and Prandtl number. To get the turbulent Reynolds numbers. Re 
— u'{L) Ljv^ we have estimated the integral scale L for each model from the wavcnumbcr fc^ for which the 
relevant energy power spectrum -E(fc) has a maximum. Together with the r.m.s. turbulent velocity and 
the kinematic viscosity, v ~ Pr(/9Cp)~^iT, we eventually obtain the Reynolds number Re and the Karlovitz 
number Ka. An example of a numerical realization is illustrated in Figure 3, where a laminar and a 
turbulent flame are shown in a space-time diagram. 

Table |l| also shows the measured turbulent flame speeds st as a function of Ka, and Pr. These 
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velocities have been obtained by measuring the slope of the function 



A{t) = / c{x,t)dx 



(13) 



where c{x, t) is the reaction-progress variable defined by (the indices b and u assign the burned and 
unburned states) 



Thus St = dA(t)/dt defines the effective flame speed. This value is constant only in the case of laminar 
flames. However, in the turbulent case it is subject to fluctuations. These estimates enable us to express the 
dimensionless turbulent flame speed, Ut = st/s\, as a function of the dimensionless turbulent fluctuation 
velocity, U = u'/si. We should mention, that we computed the dimensionless turbulent flame speed using 
the alternative formula 



estimated values are in good agreement with the ones obtained through the equation st = dA{t)/dt. 

For our studies we vary the r.m.s. turbulent velocity, the Prandtl number and the density while the 
integral length scale L remains at a constant value. 

The results arc shown in Figure 4. It clearly demonstrates that the turbulent flame speed at high 
turbulent strain rates, J7 3> 1, docs not follow a law valid in the flamelet regime. At a Karlovitz number 
greater than one, i.e. in the distributed burning regime, there is a significant deviation from this law; 
stronger turbulent intensities do not necessarily lead to a faster burning rate. This effect is more pronounced 
in the low-density case, where higher turbulent velocities even cause a slight deceleration of the flame speed. 
This result confirms the intuitive picture that fast fiames (fiames at higher density) are more resistant to 
distortion through turbulence than the slower ones (flames at lower density). 

As already mentioned in this section, the transport of heat is expected to be decoupled from turbulence 
on scales smaller than ZkPr~^^^. This assertion is numerically confirmed in the realization Bl, where we 

significantly decrease the Kolmogorov length scale (using a Prandtl number of 0.02) while keeping the 
fluctuation velocity u'(L) at a value used in the realization A7 (Pr = 0.1). The estimated turbulent flame 
speed St of both models however turns out to be nearly the same, see Figure 4. This behavior suggests 
that in cases of a small Prandtl number one should use the Peclet number, Pe = Re Pr, as the governing 
parameter for the strength of turbulent flame dynamics. 




(14) 




(15) 



where St and 5*1 denotes the turbulent and laminar profile of the energy generation rate. In fact, the 
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For clearity, Table 2 gives a summary of all the dimensionless quantities that we have used or mentioned 
in this section. 



4. Discussion 

4.1. Burning in the Distributed Regime 

The coupling of ODT as a model of homogeneous turbulence to reactive hydrodynamics gives an 
insight into the physics of distributed burning. First of all the qualitative physical picture of this burning 
process is revealed: turbulent vortices enter the reaction zone and carry away burning material to the front 
and to the back of the flame. As a consequence, one observes a completely different behavior than in the 
flamelet regime, because unlike the latter, the distributed regime changes the local flame velocity (The local 
flame velocity is the speed at which the flame propagates locally normal to itself.). 

A different feature is the existence of localized burning with high burning rates, commonly referred to 



as a local explosions ( Radford fc Chan 1995 ). Turbulent motion forms regions where hot reaction products 
are mixed with cold unburned material. These regions have a higher temperature than pure fuel and possess 
comparably short destruction timescales for carbon. Their eventual re-ignition causes high burning rates, as 
shown in Figure 5. We find that with higher Karlovitz number more of these events occur and their strength 
(burning rate) is increased. Terrestrial experiments show the generation of blast waves after the occurrence 



of local explosions. In some experiments the strongest waves developed even into a detonation (Radford & 



Chan 1995). However, laboratory experiments on turbulent flames deal with much higher expansion factors, 
Pu/pb, than those given by flames in white dwarf matter. Thus we expect the effect of blast waves in the 
context of SNe la to be less effective than in typical chemical flames. However, because our model does not 
consider pressure waves, the whole effect of local explosions on the burning process cannot be presented 
here and we cannot completely exclude the possibility that blast waves possibly enhance the global burning 
rate by compression and pre-heating of the surrounding medium. Neglecting such effects in our model 
possibly means an underestimation of turbulent combustion rate in the distributed regime. 

Taking into account the estimates on the turbulent flame speed of Table || we obtain a relatively good 
agreement with experimental data presented by Abdel-Gayed et al. (1987), henceforth AG97. Their work 
contains about 1650 experimental results on turbulent flame velocities at various strain rates U = u' / s\, 
fuels and equivalence ratios. Their results were obtained under conditions far from those in an exploding 
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Chandrasekhar-mass white dwarf, but it is worth to verify if certain physical statements, Hke the function 



Ut (Ka; Re) , appear to show universality. These experimental results and earlier work ( Abdel-Gayed & 



Bradley 1985) present two important phenomena. First, for a fixed Reynolds number, Ut is not a strong 



monotonic function of the strain rate U. Rather, it reaches a maximum at a certain turbulent intensity 
and eventually saturates or even decreases for higher values of the strain rate (bending effect) . As a second 
feature, there is flame extinction at a Karlovitz number Ka > 7 for Reynolds numbers Re > 10^. Figure 
6 shows several experimemtal values of the function Ut taken from the work of AG87. These values 
correspond to the same sample of ratios u'/si we use in our numerical caluclations. 

Our numerical results do show the bending effect for high strain rates which is in agreement with 
laboratory experiments. However, they do not show complete flame quenching, i.e. a situation where fuel 
is not consumed anymore. This is not surprising because the basic laws of thermodynamics prevent a 
turbulent (or laminar) flame from complete quenching. The first law of thermodynamics states that the 
reactive mixture must eventually reach the adiabatic flame temperature. The second law demands that 
heat must be transported from the burned to the unburned regions. Therefore, there is no way to prevent 
consumption of unburned material in a system without possible heat losses. And clearly, unlike under 
laboratory conditions, the flame in an exploding white dwarf is not subject to significant heat losses. 

The comparison with our results shows that the flame speeds presented here are lower than those in 
AG87. We think that one reason for this discrepancy is the different value of the Lewis number in our work 
and in the data of AG87. It is known that high Lewis number flames are more sensitive to strain than 
flames with a low value of Le. In fact, analytical results on planar laminar flames show that the flame speed 
decreases when a velocity field with strain, i.e. a non-vanishing gradient of this velocity field, is invoked 
parallel to the flame surface. One finds that at higher Lewis numbers the decrease of the fiame speed is 



stronger (Tromans 1981). Therefore we expect thermonuclear flames in C+0 matter to be more sensitive on 
strain effects, because of their Lewis numbers of aprrox. 10^. Another source of this deviation could be the 
already aforementioned difference in the expansion factors between laboratory and white dwarf combustion. 

Thermonuclear flames in the distributed regime thus experience a two-fold effect through turbulence: 
One is the acceleration due to an enhancement of the heat transport produced by small-scale turbulence. 
The other is a slowing-down due to additional turbulent strain, which can be measured by the Karlovitz 
number. However, we have to stress that we do not have revealed the actual physical mechanism responsible 
for the deceleration. This is the task of future work. 



-14- 



4.2. Implications for a Deflagration to Detonation Transition 

As already stated by many authors, a delayed detonation within the explosion of a Chandrasekhar-mass 
white dwarf is a reasonable model. It is a good fit to experimental constraints for element abundances 
obtained from observed spectra. In addition, although it is almost clear that after ignition the burning 
process evolves as a deflagration, there is no fundamental physical objection why a transition to detonation 
in the subsequent stages of the explosion should be impossible. Therefore, DDT models should be taken 
into account in supernova theory. In the context of this work, we discuss the possibility of distributed 
burning leading to a detonation. 

Based on the Zcldovich gradient mechanism, Khokhlov et al. (1997) calculated the minimum size, L^, 
of a region where fuel and ashes are mixed and out of which a stable detonation wave can emerge. In their 
model the formation of this region is caused by turbulent mixing, but it is not obvious how this process is 
realized. They proposed that global expansion of the star along with the following density drop quenches 
the deflagration flame, whereas the Rayleigh- Taylor instability still takes place on scales of 10^-10^ cm. 
Before the star recompresses again, fuel and ashes are mixed and turbulence eventually produces a complete 
mixed region of size of the order of Lc- Thus, pulsation may lead to a DDT. 

Apart from considering pulsations of the star, it is worth asking whether burning in the distributed 
regime alone can lead to a DDT. Niemeyer & Woosley (1997) advocated this burning mode as a favorite 
model for DDT. Basically, they proposed that if turbulent burning is able to form a region of uncompletely 
burned material with a certain temperature gradient, then a carbon detonation can occur when two 
presumtions are fulfilled. First, the aforementioned region must have at least a size of the critical radius 
necessary for a detonation (These critical radii could be calculated explicitly for different densities and 
nuclear compositions, see Niemeyer & Woosley 1997.). Second, the temperature gradient must be shallow 
enough in order to allow a sonic phase velocity of nuclear burning. For instance, consider a mixed region 
consisting of half carbon and oxygen at a density of 3 x lO^g cm~^ and having a maximum temperature of 
T'max = 2 X 10® K. It then follows that this region must have a size of at least Lc ^ 50 m. Furthermore, 
temperature differences within it should nott exceed 10^ K. Then a stable detonation wave can form. 

Both approaches have in common that they require a rather shallow temperature gradient over a region 
with a sufficient amount of carbon. Until now, it has not been demonstrated how turbulence can achieve or 
lead to such a high level of temperature homogeneity. 

In this context it is interesting to know the maximum length scale, Zmax, where distributed burning in 
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a white dwarf can occur. A lower bound of this scale is given by the size of those eddies whose turnover 
time equals a typical turbulent burning timescale Tnuc- That is 

^max/^ (^inax) — '^nuc ■ (-^^) 

In the case of laminar flames Tnuc is just the diffusion timescale t\ « S\/s-i « 0.7 x 10^^ s. But what is 
the value Tnuc in the presence of turbulence? We argue that even in the turbulent case the choice of ti as 
the typical burning timescale is natural. To justify this assumption, we compute the spatial distribution of 
destruction timescales across the flame. In general, this distribution is characterized by a minimum value 
T*n^. It corresponds to the location of the maximum energy generation within the flame. In our numerical 
models, we find that t*^^ is independent of the turbulent intensity, i.e. its value is the same for the laminar 
case and for the distributed regime. In fact, t*^^ equals t\. Then equation (|l^) leads to the useful expression 

'max~^lKa. (17) 

This number can be viewed as the size of the distributed flame, i.e. an estimate on the size of the interface 
that seperates burned and unburned material. However, there are reasons to believe that in reality Z^ax 
exceeds the value given by equation ([l7|). This is because there are turbulent vortices with bigger size than 
Zmax that can enter the flame and transport burning material with longer burning times than Tnuc- Hence, 
we do justice to these considerations and set the effective value of Zmax to 

'max ~ a Ka , (18) 

with a > 1. The actual value of a however, i.e. the ability of turbulence to stretch out the flame brush, has 
to be estimated by more rigorous and quantitive considerations, see Lisewski et al. (1999). 

If 'max itself reached the critical size Lc required for a detonation, then burning in the distributed 
regime would become more of a promising model for DDT. But still it would not be clear if present 
temperature fluctuations would be small enough to allow this transition to happen. Thus the equality of 
'max and Lc is only a necessary condition for a DDT. At a density of 2.3 x lO^g cm^"^, Khokhlov et al. 
(1997) obtain Lc ~ 10'' cm. Now the expected turbulent velocities at an integral scale of L = 10^ cm are 
u' ^ y/9cffL ^ lO^cm s^^, where gcS ~ lO^cm s^^ is the effective gravitational acceleration. Hence, at 
densities p ^ 2 x lO^g cm^"^ the corresponding Karlovitz numbers are only around Ka « 10. According 
to equation ( p^ a very high Karlovitz number Ka w 10"* is required to obtain equality of 'max and Lc- 
Consequently, the enhancement-factor a in equation (|l^) would have to be of the order of 10'^. This means 
that the ashes- fuel interface has to be extended much more in size than given by equation (p^). 
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Also at a density of 10''g cm where s\ = 4.7 x lO'^cm s ^ and S\ = 4.2 cm (Timmes & Woosley 1992), 
the value of a has still to be around 10"^ in order to enable a direct DDT. 

Furthermore, the process of a DDT demands, that the mixing timescale, r^i^, which can be regarded 
as the time it takes turbulence to form a complete mixed region of size Lc, must be smaller than a typical 
burning timescale, t^uc, of the mixed turbulent matter. But in order to avoid quenching by expansion, t^^c 
in turn must be smaller than the hydrodynamical time of the star, Th, c.f. equation (^). These conditions 
read as Tmix < Tnuc < Th, or as 

Lc 446 
u'(Lc) p^l^ 

For p = 2 X lO^g cm~^, we have Tmix ~ 5 ms while the hydrodynamic time is Th = 0.02 s. Here now the 
question arises whether turbulence can form such a region whose burning time is smaller than r^. This 
problem is adressed to forthcoming work, where we study - by using ODT - the burning times on different 
lengthscales within the distributed flame. However, until we do not have a better understanding of the 
large scale structure (Icnght-scales of ^ Lc) of flames in the distributed regime, the preliminary result from 
equation suggests that Karlovitz numbers in white dwarf C+0 matter at densities around lO^g cni^'^ 
are too small to form a region of critical size necessary for a DDT. 



5. Conclusions 

The aim of this work was to investigate the physics of distributed burning in an exploding white dwarf 
of chandrasekhar mass. With ODT as a new model for homegeneous, isotropic turbulence we are able to 
study the effective turbulent flame speeds resulting from this kind of burning in dense C+0 matter. Our 
results show that in this regime the local properties of the thermonuclear flame are readily changed: the 
flame cannot be represented by the flamelet model anymore and moreover, higher turbulent intensities in 
general can lead to a lower local flame speed. This behaviour has already been observed in laboratory 
experiments and is introduced the first time in the context of nuclear explosions in type la SNe. We have 
also found that, in cases where the Gibson length Iq, is smaller than 8\ and where the Prandtl number 
is much smaller than one, turbulent burning decouples from turbulent scales smaller than /kPr^'^^^. This 
suggests that the appropriate parameter for turbulent heat transport inside a white dwarf is given by the 
Peclet number instead of the Reynolds number. Because of the small ratio of viscous momentum to diffusive 
heat transport the former is by several orders of magnitude smaller than the Reynolds number: Pe = Pr Re. 

In addition, first quantitative considerations demand that a turbulence-induced deflagration to 
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detonation transition can only happen deep in the distributed burning regime, i.e. at Karlovitz numbers 
larger than 10^. However, our estimates for the turbulent velocity at the integral scale in an exploding 
white dwarf give values of Ka that are too small (at least by a factor of 10^) in order to make a direct DDT 
probable. 

Finally, it is wortwhile to remark that we have completely left out questions about the influence of the 
turbulent spectrum on this issue. We have only shown that a DDT is unlikely, when Kolmogorov theory of 
turbulence is considered. Niemeyer & Kerstein (1997) proposed that turbulence in an exploding white dwarf 
is dominated by a potential energy cascade instead of Kolmogorov's kinetic energy cascade. This domination 
leads to different scaling properties (Bolgiano-Obhukov scaling) of the turbulent velocity. As a consequence, 
higher turbulent velocities are necessary for the transition to distributed burning. Furthermore, Niemeyer 
& Woosley give the speed of the turbulent fluctuations, urt ~ O-^VgesL « u'{L), which may not be the 
upper limit for turbulent velocities at the largest turbulent scale. It is based on the original idea of Kerstein 
(1996), that the turbulent flame velocity, st is increased by local expansion caused by thermonuclear 
reactions {active turbulent combustion). At the largest length scale, L, and at a certain temporal stage of 
the burning process, the speed of the flame may exceed urt- Because the relation u'{L) w st still holds at 
this point, a significant enhancement of turbulent velocity fluctuations may occur. But it is unclear, if the 
growth rate, st/st, is high enough to realize this enhancement within a hydrodynamical time Th. 
These arguments favour a DDT, because their effect is an increase of the turbulent fluctuations on lengths 
comparable to L,,. However, both arguments arc rather speculative, because they have not been confirmed 
numerically, analytically or experimentally in the context of SNe la yet. 
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Table 1. Realized models of turbulent flames. X(^^C) = X{^^0) = V2, p = 2.3 x 10^ g cm'^ (Al - A9, 

Bl), p = 1.3 X lO^g cm-3 (CI - C5, C41, C51). 



Model 


Reynolds Number 


Karlovitz Number 


Prandtl Number 


Flame Speed 




Re 


Ka 


Pr 


St fxlO'* cm s^^l 

^ L J 


Al 


114 


0.06 


0.1 


2.15 


A2 


270 


0.23 


0.1 


2.4 


A3 


452 


0.51 


0.1 


2.7 


A4 


598 


0.78 


0.1 


3.1 


A5 


673 


0.92 


0.1 


3.4 


A6 


820 


1.24 


0.1 


4.2 


A7 


1285 


2.43 


0.1 


5.0 


A8 


1493 


3.05 


0.1 


5.3 


Bl 


6182 


2.29 


0.02 


4.9 


CI 


28 


0.68 


0.1 


1.0 


C2 


119 


5.80 


0.1 


1.8 


C3 


212 


13.9 


0.1 


2.4 


C4 


314 


25.7 


0.1 


2.1 


C41 


323 


25.3 


0.1 


2.2 


C5 


366 


32.4 


0.1 


1.9 


C51 


371 


32.2 


0.1 


2.0 
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Table 2. Summary of some important dimensionless parameters. 



Symbol 



Name 



Mathematical definition 



Physical description 



Re Reynolds number 



Pr Prandtl number 



Pe Peclet number 



Ka Karlovitz number 



Le Lewis number 



u'{L) L 

V 

a 

Re Pr 

Silk 



pCpD„ 



integral scale turbulent difFusivity 
viscous diffusivity 

viscous diffusivity 
heat diffusivity 

integral scale turbulent diffusivity 
heat diffusivity 

heat diffusion timescale 
smallest turbulent heat transport timescale 

heat difFusivity ^ 
mass diffusivity 



^Dui denotes the ass diffusion coefficient. 
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0.1 1.0 10.0 100.0 

k [cm-^] 

Fig. 1. — An energy power spcetrum generated by stationary ODT (o). Up to a constant factor the power 
spectrum is the Fourier transform of the second structure function 5*2 — {{u{r) — u{0))^). The sohd hne 
represents the k~^^^ law for the inertial range predicted by Kolmogorov theory. 
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Fig. 2. The propagation of a X(^^C) = X(^^O) = 0.5 conductive flame. The sequential temperature profiles 
(dotted lines) are shown at times t = 1.3 x lO""*, 3.0 x lO"'^, 6.5 x lO""*, 8.3 x lO"'', 1.01 x IQ-^, 1.19 x 
10-3, 1.36 X 10-3, 1.51 X 10-3 s. The dashed line is the density profile at t = 1.19 x 10-^ s. The carbon 
flame already moves in a steady-state, while the temperature rise at the left-hand-side comes from oxygen 
burning due to the energy stored at the initial state of the simulation. 
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Fig. 3. — Comparison of a laminar flame (left panel) propagating into C+0 matter with a turbulent one 
(right panel) generated with ODT. The lines represent iso-temperature regions at values from 0.6 to 3.0 x 10^ 
K. The right panel shows data from the numerical realization A2. 
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op = 2.3x lO^gcm^ 
A /) = 1 .3 X 1 0^ g cm~^ 
XC^C) = XC'O) = 1 /2 
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Fig. 4. — The dimensionless turbulent flame speed as a function of the dimensionless r.m.s. turbulent 
velocity fluctuations. At high strain rates there is a strong deviation from the flamelet model, which predicts 
Ut — {1 + (3U^y/^ (solid curve with (3 — 0.24), (see for instance Pocheau 1994). The numerical realization 



Bl is depicted by a box (□). The two black triangles give the results from the realizations C41 and C51 
having a four-times higher spatial resolution than C4 and C5. The gray triangles represent those speeds 
obtained from formula (|l5|). 
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Fig. 5. — A local explosion at the front of a turbulent flame in the distributed regime. The profiles are taken 
from the model A3 
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Fig. 6. — Experimental data for the turbulent flame speeds taken from Abdel-Gayed et al. 1987. 



